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CORRELATORS AND CROSS-CORRELATORS 
USING TAPPED OPTICAL FIBERS 



REFERENCE TO RELATED APPLICATION 
This application is a continuation-in-part of U.S. Patent Application Serial No. 
09/765,215, filed January 18, 2001, the entire content of which is incorporated herein by 
reference. 



5 FIELD OF THE INVENTION 

This invention relates generally to signal processing and, in particular, to signal 
correlation and cross-correlation using tapped optical fibers. 



BACKGROUND OF THE INVENTION 
As fiber optic technology continues to grow, the use of various diffraction 

10 gratings within the fibers for signal processing will become increasingly common. Our 
co-pending U.S. Patent Application Serial No. 09/765,215, incorporated herein by 
reference, describes a two-dimensional, fiber-delay radiator (FDR). The device is made 
with optical-fiber taps forming Bragg gratings orientated at 45 degrees to the fiber core, 
permitting light to emerge directly out of the side of the fiber, linearly polarized. The 

15 optical fiber is wound on a cylindrical-like form such that a number of loops of the fiber 
are available for making a number of taps on each loop. Taps are preferably generated 
along each loop of the fiber so that a portion of the light propagating in the fiber will exit 
sideways from fiber at the taps. A lens system is then used to capture the light from the 
taps and produce a Fourier Transform of the total distribution of light from all the taps. A 

20 video camera then captures this Fourier Transform light and the video signal is acquired 
for further computer processing. The result is that the power spectrum of the light signal 
is displayed on a monitor. 

The preferred construction of the FDR would include taps with "ideal" phase 
characteristics. It may be possible to directly fabricate such a device. However, a 
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method to correct the phases in real-time and to introduce desired weighting may be 
preferable. To achieve this goal, a phase spatial light modulator may be used to correct 
and modify the tap phases. Another embodiment uses a coherent reference wave to 
generate a holographic optical element (or complex spatial light modulator) to correct the 
5 tap phases and amplitudes. According to this technique, a coherent reference source and 
a detector array are used to capture the radiation amplitude pattern of the FDR. Then, 
with digital processing of the captured pattern, the desired spectral signal properties are 
obtained. 

This basic configuration finds application in a wide variety of commercial, 
10 industrial and military situations including, but not limited to the analysis of signals that 
have been modulated on light propagating in a fiber; the analysis of telecommunication 
channels in optical fiber network links; the analysis of optical sources by determining 
their spectral characteristics; and as part of a spectrometer to analyze gases and other 
materials. The spatial modulator version, in particular, could be used as a phased array 
15 for Ladar and Lidar applications; as a beam director for data storage and other scanning 
applications; and as a channelizer for use in separating out and monitoring individual 
signal channels. 

SUMMARY OF THE INVENTION 
This invention extends the principles disclosed in co-pending U.S. Patent 
20 Application Serial No. 09/765,215, with emphasis on the use of a Tapped Optical-Fibers 
Processor (TOP) for correlating wide bandwidth and large time-bandwidth signals. In 
particular, it can be applied to processing radar and SAR (synthetic aperture radar) 
signals. The TOP can also be used for lidar and ladar systems with out the need for 
modulators or separate carrier light source. Broadly, using this approach, fine resolution 
25 can be obtained without fast front-end sampling and digital computational burdens can be 
significantly reduced. 

A preferred method of signal processing according to the invention includes the steps 
of providing an optical carrier signal; modulating the optical carrier signal with an input signal 
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to provide an optically modulated signal; radiating the optically modulated signal from a set of 
taps formed in an optical fiber; performing a spatial Fourier transformation on the radiated 
signal; detecting the Fourier transformed signal and converting the detected signal into an 
electrical signal; and performing a digital Fourier transformation as well as other processing 

5 functions on the electrical signal to output an autocorrelation of the input signal. 

The radiated signal is a function of the distance between the taps, the velocity of the 
signal through the fiber, and an aperture weighting function. The resulting the autocorrelation 
of the input signal is weighted by the autocorrelation of the aperture weighting function. 

Particularly in conjunction with radar signal processing, the input signal may be 

10 composed of the sum of at least two or more signals, in which case the output may include the 
autocorrelations of both inputs as well as the generation of cross-correlation of the two signals. 

In terms of hardware, a signal processor according to the invention preferably includes 
a coherent laser source operating at a carrier frequency; a modulator to insert an input RF 
signal into the carrier; an optical fiber radiator composed of a fiber with taps that radiate the 

15 modulated optical signal; a lens to perform a spatial Fourier transformation on the radiated 
signal; and a detector array to output the transformed signal to a digital processor for 
additional signal processing. In any case, the two input signals may be electronically or 
optically combined. 

In the case of ladar and lidar applications, the received light along with the source light 
20 could be used as input to the invention without the need for modulators or a separate carrier 
light source. 



BRIEF DESCRIPTION OF THE DRAWINGS 
FIGURE 1 illustrates a Tapped Optical-Fibers Processor (TOP) system applicable 
to the present invention; 
25 FIGURE 2 shows how the tapped Optical-Fiber Processor (TOP) can be used in a 

ladar or lidar system; 

FIGURE 3 shows the correlation outputs for dual signal inputs with delay offsets; 
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FIGURE 4 depicts front-end usage of the TOP as part of a synthetic aperture radar 
system using a spotlight processing algorithm along with alternative processing 
techniques; 

FIGURE 5 shows a TOP configuration with dual modulators and variable fiber 
5 delay; and 

FIGURE 6 presents a bandwidth versus aperture delay comparison between 
commercially available AO cells (Xs and Os) and potential FDRs (diamonds). 

DETAILED DESCRIPTION OF THE INVENTION 
The basic TOP system concept is illustrated in Figure 1. It consists of a coherent laser 
10 source operating at the optical carrier frequency v c , a modulator to insert RF signals onto the 

optical carrier, an Optical Fiber Radiator (OFR) composed of a fiber with taps that radiate the 
optical signal, a lens to take the spatial Fourier Transform of the radiated signal, and a detector 
array that read out the transformed signal to a digital processor for final processing. 

First consider a RF signal u(t) which is inserted onto the optical carrier signal to give 
1 5 the modulated optical signal 

s(t) = u(t)e- 2 ^ (1) 
This signal is radiated from the taps in the FDR as it propagates along the fiber. The radiated 

light is given by 

a(x,t) = s(t-p-w(x) (2) 

20 where x is distance along the FDR taps, V is the effective velocity of the signal along the taps, 
andw(;t) is the aperture weighting function. This radiated signal is then spatial Fourier 
Transformed by the lens to give the following optical amplitude signal at the detector array 

A(fj) = S(f)e- 2 " itf *W(f) 

= U{f)e- 2nitf *W{f)*S(f-v c ) 

where / is a normalized distance with units of frequency, * represents the convolution 

25 operation and upper case letter represent the forward FT function of functions represented by 

4 
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corresponding lower case letters. (Note: / = - — , where L is the focal length of the lens, X c is 

the operational optical wavelength, and z is the physical coordinate with the same positive 
sense as for x.) 

The TOP is set to operate with the optical carrier at the origin where / = 0 . The 
5 instantaneous intensity that is seen by the detector array is thus given by 

I(f,t) = \A(f-v c ,tf 

(4) 

= \U(f)e-™' f *W(f)\ 2 
The detector array then integrates the intensity to generate the acquired signal 

t n +T d /2 f 9 +Ff2 

K q =V d J dt \ dfl(fj) q = l,...,Q (5) 

t n -T d /2 f q -F/2 

where I is the signal for qth detector at acquisition time t n , T d is the time integration, F the 
10 spatial integration, and rj d the conversion efficiency. 

We assume that the acquisition sampling in space is adequate to represent the 
intensity function as a continuous function. This allows us to use Fourier Transforms rather 
than discrete digital Fourier Transforms in the following analysis. Thus, the detected signal 
can be written as 

15 /„(/) = rj d ]rect(^)I(fj)dt (6) 

and its FT as 

K W = V d \ rect£yiL)i(T,t)dt (7) 

— ee (j 

where the FT of the intensity signal I(f,t) is 

i(r,t) = [u(r-t)w(T)]'®[u(r-t)w(r)] . (8) 

T 

20 Here, ® represents the correlation operation in delay space as defined by 
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Thus, in general, the correlation is between weighted traveling waves of the input signal. 

Now consider the detected signal captured by the array by substituting Eq. (4) into Eq. 
(6), which gives 

5 /„(/)=% ldf'\dfv\f')u(rw\f-f'w(f-rf^ do) 

Carrying out the time integration gives 

'.</> = Vj d \df \dfV*(f')U(nW* (f - f')W(f - /") sinc[r, . (/' - f)¥** um - f ** 

(11) 

We now note that if our integration is long with respect to any delays of interest, the 
10 normalized sine function can be replaced by a delta function and we can let /' = f* to obtain 

I n (f) = V d \df\U(f)f\W(f-f)\ 2 

= rj d \U(f)\ 2 *\W(f)f 
We now take the FT with respect to /to obtain 

i n (T) = TJ d [u(T)®u(TMw\T)®w(T)] (13) 

This is our basic result, namely, that the FT of the detected signed gives the gives the 
1 5 autocorrelation of the input signal weighted by the autocorrelation of the window Junction. 
If the aperture function is a uniform weight given by 

w(r) = rect(— ) , (14) 

where T a is the transit aperture time, then 

w*(r)®w(r) = A(— ) (15) 

20 where the triangle function is defined by 

l-|r| if |r|<l 



A(r) = 
6 



0 if|r|>l 



(16) 
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We note that T a = X a I V , where X a is the aperture length of the FDR. It is clear, physically 
and mathematically, for a continuous signal that as the delay r increases the common overlap 
in signal decreases and no correlation is obtained for |r| > T a . 

Cross-correlation Operation 

5 Of particular interest is when the input RF signal is composed of the sum of two 

signals. For this case of interest the input signal is given by 

u(t) = u a (t-T a ) + u b (t-T b ) (17) 

Then, the signal correlation from Eq. (13) is given by 

U{T)®U(T) = C{T) (lg) 
= C aa (T) + C bb (T) + Crt(T-AT) + C ba (T + AT) 

10 where the correlation functions are defined by 

c m (t) = u m {T)®u n {T) m,n = a,b (19) 
and Ar = r b -r a . Thus, as indicated in Figure 3, the cross-correlations are offset by their 
delay differences from the centered auto-correlations. 

We now assume, as is typical for radar or ladar or lidar signals, that if u a (t) = u 0 (t) is 
15 the transmitted signal then the received signal can be represented by 

u b (t) = u 0 (t)f d b m S(t-T m ) (20) 

m=l 

where b m is the radar amplitude from a point scatter with relative time delay r m , and M is the 
number of scatterers. Or we can describe the return in terms of the scattering distribution 
b(t) as 

20 u b (t) = u 0 (t)®b(t) (21) 

Thus, we can write 

c(r) = C(m it) + c bb (r) + c ab (r - A r) + c ba (r+ At) (22) 
where for point scattering description 
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c flfl (T) = c 00 (r) = « 0 *(r)®« 0 (r) (22) 



M M 



c bb (T) = YZKb n 5(j-T m +T n ) (23) 

m /i 

M 

^W = c 00 (T)£fo>(r-rJ (24) 

m=l 



m=l 

5 or for a scattering function description 

= kw® ® k W® W] (26) 

^(r) = c 00 (r)O6(r) (27) 
c ta (r) = 4(T)®6(r) (28) 
It is seen that the scattering object is range imaged in the cross-correlation terms with 
10 the autocorrelation term c m (f) as the point spread function. See Figure 3. If the signal has a 
bandwidth of 5, then the delay resolution is 

If &T range is the maximum spread of delays in the range signal u b (/) , then we see that 

the following condition needs to be satisfied for the image to be captured without exceeding 
15 the aperture window or overlapping the autocorrelation terms. 

3Ar Ar 

SL<Ar<r =s- (30) 

2 2 

The maximum delay spread is given by 

MoxA W =|- (31) 

when At = 37;/ 2. 
20 Configurations 
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The TOP system would be used as indicated in Figure 4. It would be part of the front 
end of a radar receiver that uses two antennas. The output of the TOP would then be further 
digitally processed. The signals from the antennas could be combined electronically to 
generate the input into TOP system, as shown in Figure 1. Alternatively, the separate signals 
5 could be combined optically as indicated in Figure 5. Here, a fiber coupler is used to separate 
the input light into two beams that are individually modulated. Then one of the modulated 
outputs is delayed appropriately via a variable fiber delay line to bring the two signals within 
the delay aperture of the processor and with delay differences to optimize the processing 
aperture. The fiber delay line could be a combination of discrete long delays (microseconds) 

10 for any fixed offsets and fine control to change the delay as a function of time, as is useful for 
synthetic aperture imaging. The two modulated signal are then recombined via a fiber coupler 
to feed the FDR. For spotlight processing the signal from the TOP could be used directly, after 
digitization, as phase history input. For back-projection processing, the digitized signal would 
be fast Fourier transformed first to bring it into the correlation (projection) domain, before 

1 5 further processing. Direct optical processing is also possible. 

It is seen that IF frequency mixing and fast digitization of the input signal is avoided 
by applying direct RF modulation of the light, which then produces a slowly changing 
interference pattern. Pre-summing of the signal is accomplished via the detector integration of 
the interference pattern. The interference pattern can be digitized at moderate rates. In other 

20 words, the TOP system has analog processed the signals, considerably reducing the required 
(subsequent post-optical) sampling rates. 

Potential Performance and Available Hardware 

Depending on physical size limitations, a range of processor parameters are available. 
Consider the following system that probably is near the limits of current technology. The 
25 primary fixed parameter is the velocity of light (fundamental mode) in a fiber, which is about 
2xl0 8 m/sec, or about 5 nsec delay per meter of fiber. If we use a 125 fim diameter fiber, we 
can get about 7 loops per mm of FDR height, and with a 20 cm high FDR, we get about 1400 
loops per FDR. If we set the loop length to 10 cm, we get a sample rate of 0.5 nsec, which 
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gives a device bandwidth of 1 GHz. Modulators are available with over 40 GHz of bandwidth, 
and hence are not the limiting factor. The total fiber length is 140 m, which gives an aperture 
time of 700 nsec. Since the FDR could be less than 2 cm wide, we could stack 10 of them in a 
20 cm wide package to get a total aperture time of 7 fisec. 

5 The resolution of the system is given by the autocorrelation of the signal. For a 1 GHz 

bandwidth signal, the delay resolution would be 1 nsec, which corresponds to half foot 
resolution for a monostatic system. 

The detector array pitch must be sufficient to sample the highest fringe frequency at 
greater than twice per cycle. The total number of detectors must be equal to twice the desired 

10 number of scene resolution elements. Of course, some overhead in sampling and number of 
elements is desirable. We typically run the detector operation near saturation to optimize the 
signal-to-noise ratio and minimize the temporal sampling rate. This is done by controlling the 
integration time and light level. The detector array frame rate is independent of the input 
bandwidth. However, for synthetic aperture radar (SAR) operation the detector integration 

15 time is determined by the azimuth sampling requirements, which tend to be a few thousand 
hertz. Hence, we might nominally require a 10 KHz frame time for the detectors, which we 
note is easily accomplished for linear arrays. Linear detector arrays are available with 44 KHz 
frame times for 4096 element arrays and correspondingly faster frame rates for fewer 
elements. In this case, the overall light level will be increased to operate near saturation. 

20 We have demonstrated basic spectral analyzer operation down to 10 MHz resolution 

and up to 100 GHz bandwidth operation at a wavelength 840 nm. Modulated spectral analysis 
up to 10 GHz was demonstrated. This spectral analysis work has used holographic techniques 
to generate digital matched filter correction of the taps. Hence, the no attempt was made to 
generate diffraction limited optical quality of the tap array, as would be desirable for this radar 

25 correlator operation. We expect that such digital filtering or an electronically controlled spatial 
light modulator (e.g., a liquid crystal spatial light modulator) will be incorporated as part of the 
FDR to give the desired tap phase correction and allow for adaptive adjustments in the analog 
optical processing. 

10 
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Acousto-optic (AO) cells have been used for correlation processing for many years. 
Figure 6 shows a comparison in bandwidth and delay aperture between AO cell technology 
and potential FDR technology. The Xs and Os represent commercially available AO cells. 
The diamonds represent FDR potential. The left two diamonds represent capabilities that can 
5 be obtained with our current components. The upper right two diamonds represent potential 
capabilities that appear feasible. It is seen that the much higher bandwidth capability exists for 
the FDR devices, even at microsecond time delays, then from AO cells. Also the FDR does 
not require the bulk optical illumination hardware, such as lenses, that is required for AO cells. 

Technology Development 

10 There are a number of aspects of the TOP system that require further analysis and 

development. 

One primary issue is the quality of the tap array. Ideally, we would like to generate a 
high optical quality tap array with prescribed positions, efficiencies and phases. The technique 
we have developed for generating the taps allows tap efficiency to be monitored while a tap is 

15 being created and to stop the tap generation at the correct level. However, the current FDR 
devices have been made with standard fibers that do not maintain the polarization state. In 
these fibers, the polarization state of the light varies as it propagates along the fiber. The taps 
produce linear polarization light, but its intensity will vary depending on the local state of 
polarization. This polarization state depends on the frequency of the light and temperature of 

20 the fiber. Using polarization maintaining (PM) fibers that will allow the generation of taps 
that will maintain their output intensities. Using the PM fiber leads to a better overall FDR 
efficiency since we are tapping a single polarization state and not losing energy to the 
orthogonal state. 

Obtaining taps with phase fidelity to a fraction of the optical wavelength is difficult. 
25 Further, maintaining the tap phases of the array will be required. The TOP center frequency 
can drift with temperature changes. Construction techniques can be used that compensate for 
temperature effects, as well as controlling the temperature of the TOP. However, two-phase 
compensation techniques are available for obtaining phase fidelity without having geometric 

11 
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fidelity. One technique is to use digital holographic recording techniques and apply a digital 
matched filter in the post processing to correct for the original fabrication errors as well as 
drifting. This holographic technique has been demonstrated. The other technique uses a 
spatial light modulator (SUM) under electronic feedback control to continually correct the 

5 output phase of each individual tap. The SLM can compensate for the original fabrication 
eirors as well as drifting. Further, it can apply real-time adaptive techniques on the array to 
compensate for dynamic effects such as could be generated by atmospheric effects. The 
holographic technique would require additional post-optical digital processing. The SLM 
technique would cause more light inefficiencies and introduce additional analog hardware 

10 complexities. 

The application of a lenslet array to the tap array would improve overall light 
efficiency by one or two orders of magnitude. It also could improve the compactness of the 
correlator system. 

Using the TOP at the front end of a receiver where the dynamic range requirements 
15 are small helps avoid pressing on the typical limited dynamic range of an optical analog 
processing system. 

Conclusions 

The TOP (Tapped Optical-Fiber Processor) clearly has unique potential to improve 
throughput rates and/or packaging when compared to an all-digital package. Extensive 
20 continuing investment by the optical telecommunication industry into fibers, fiber based 
devices and techniques, generates an improving technological base for successful 
development of this new type of fiber-based correlator. 
I claim: 
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